Metal-insulator transitions in the Periodic Anderson Model 
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We solve the Periodic Anderson model in the Mott-Hubbard regime, using Dynamical Mean Field 
Theory. Upon electron doping of the Mott insulator, a metal-insulator transition occurs which is 
qualitatively similar to that of the single band Hubbard model, namely with a divergent effective 
mass and a first order character at finite temperatures. Surprisingly, upon hole doping, the metal- 
insulator transition is not first order and does not show a divergent mass. Thus, the transition 
scenario of the single band Hubbard model is not generic for the Periodic Anderson model, even in 
the Mott-Hubbard regime. 
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The metal-insulator transition in strongly correlated 
materials remains a central problem of modern condensed 
matter physics 0, 0] • Great progress in its understand- 
ing was made possible by the development of new the- 
oretical approaches such as the Dynamical Mean Field 
Theory [3[ , which is a method that becomes exact in the 
limit of large lattice connectivity The mean field 
equations can usually be tackled with a variety of nu- 
merical approaches which allow to obtain reliable solu- 
tions and insights. In this context, the Hubbard model, 
which is probably the simplest model that captures a cor- 
relation driven metal- insulator transition (MIT), called 
Mott-Hubbard transition, has received most of the at- 
tention in the past 15 years. As a result of intense inves- 
tigation, our understanding of the metal-insulator tran- 
sition in that model is now profound. The studies have 
unveiled a scenario where, at low temperatures and mod- 
erate interaction, the half-filled Mott insulator may be 
driven to a correlated metallic state through a first order 
transition [Bj]. The transition can occur as a function of 
correlation strength, temperature or doping. The first 
order line ends at finite temperature in a critical point 
and the critical region can be described by a Ginzburg- 
Landau theory This theoretical prediction was ex- 
perimentally verified in experiments on V2O3 [?]. The 
Hubbard model is often considered as a minimal model 
for the study of rather complicated compounds such as 
transition metal oxides and heavy fermion systems. This 
is supported by the implicit assumption that the Hub- 
bard model is expected to be the effective low energy 
Hamiltonian for a wider class of more realistic multiband 
models for strongly correlated electron systems. 

On the other hand, a more realistic model which is 
also widely used in theoretical investigations of strongly 
interacting systems, though still schematic, is the Peri- 
odic Anderson model (PAM) . In the context of correlated 
electron systems, this model permits to describe explic- 
itly both, the localized orbitals, such as the d in transi- 
tion metal oxides or the / in heavy fermion systems, and 



their hybridization to an itinerant electron band (such 
as that of p orbitals of oxygen in transition metal ox- 
ides). In fact, the PAM allows to investigate the various 
regimes where Mott insulating states occur, as charac- 
terized by the Zaanen-Sawatzky- Allen (ZSA) scheme [||. 
They are classified as either Mott-Hubbard insulators or 
charge transfer insulators. The first apply to the early 
transition metal oxides such as titanates and vanadates, 
while the second is relevant for cuprates, such as the high 
T c superconductors,and manganites, which show colossal 
magnetoresistance [2|. In theoretical studies, however, it 
is often assumed that both Mott-Hubbard and charge 
transfer systems may be described at low energies by a 
simpler one band Hubbard model Hamiltonian. 

In the present work we shall test the putative validity 
of the Hubbard model as the effective low energy Hamil- 
tonian of the more realistic Periodic Anderson model. 
We shall do this within a well defined mathematical 
framework, namely, the Dynamical Mean Field Theory 
(DMFT), that allows us to obtain essentially exact nu- 
merical solutions of the models (in the statistical Monte 
Carlo sense). In particular we shall concentrate on the 
nature of the (paramagnetic) metal-insulator transitions 
that occur in the Periodic Anderson model with param- 
eters that set it in the Mott-Hubbard regime, and dis- 
cuss it with respect to the corresponding scenario that is 
realized in the one band Hubbard model case. In addi- 
tion, our results should also be valuable for the interpre- 
tation of experimental spectroscopies of strongly corre- 
lated transition metal oxides, that experienced fantastic 
improvements in the last decade. In fact, the analysis of 
experimental data of systems which have a mixed orbital 
character is not always simple when strong correlations 
arc present. Finally, our work addresses a very relevant 
issue in regard of the intense effort that is currently ded- 
icated to the implementation of ab initio methods for 
strongly correlated materials M which makes heavy use 
of the DMFT methodology [lj. 

Among our main results we find that in the case of the 
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electron doped driven MIT, the scenario is indeed simi- 
lar to the one realized in the Hubbard model, however, 
the hole doping scenario is qualitatively different. In this 
case, the correlated metal has a resonance peak at the 
Fermi energy flanked by a Hubbard band, but, unlike 
the Hubbard model scenario, it is not related to the for- 
mation of a Kondo like resonance and its mass does not 
diverge at the transition. Moreover, and also in contrast 
to the Hubbard model case, our results indicate that this 
metal-insulator transition is of second order as no signs of 
coexistent solutions were observed. We shall argue that 
while the metallic state in the former case is a renormal- 
ized "Brinkman-Rice" Fermi liquid the latter can be 
interpreted as liquid of "Zhang-Rice singlets" [111 ]. 
The Periodic Anderson model Hamiltonian reads, 



cavity method 0, , and one obtains the local effective 
action for the electrons: 
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where the d a and operators destroy and create elec- 
trons at non-dispersive d— orbitals with energy ed, p a and 
p+ destroy and create electrons at p— orbitals with energy 
e p which form a band with hopping parameter t. The p— 
and d— orbitals are hybridized with an amplitude t pd , and 
the electron correlations are introduced by the Coulomb 
interaction U on the d— sites. It is customary to define 
the charge transfer energy A = e d — e p , and fi is the 
chemical potential. As described in the ZSA scheme this 
model predicts correlated insulating states in two very 
different regimes: at A << U the charge transfer insu- 
lator, and at U < A the Mott-Hubbard insulator. The 
latter is relevant for the early transition metal oxides and 
will be the focus of the present work. 

To solve the PAM using DMFT, for simplicity we adopt 
a Bcthe lattice that corresponds to a semicircular density 
of states (DOS) for the p— electron band. Setting the hop- 
ping of the p— electrons to t — 1/2, their half-bandwidth 
is equal to one, and fixes the units of the model. To 
set the system in the Mott-Hubbard regime, we adopt 
€d = and e p negative, so that the p— band lies well 
below the Fermi surface and is almost full, while the oc- 
cupation of the local d— sites will be close to one. The 
parameter t p d controls the hybridization between the or- 
bitals at each lattice site, and permits the derealization 
of the d— electrons. In fact, a finite t p d turns the "flat" 
band of d— orbitals into a conduction band with mainly 
d character and bandwidth of the order of t pd /A. Now, 
for a moderate value of the repulsion U > tp d /A and 
an occupation of the d— site n d close to one, one expects 
the conduction band to open a correlation gap and the 
system becomes a Mott-Hubbard insulator. 

The DMFT equations are most easily derived using the 
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where d a and d+ correspond to a given (any) site of the 
lattice. This equation defines the so called associated 
impurity problem of the model, that is subject to a self 
consistent constraint that reads, 
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The solution of the quantum impurity problem |T]) gives 
the local d— electron Green's function G dd and defines a 
self-energy Y,=Q^ 1 — G dd . The local Green's function of 
the p— electrons G pp is obtained in terms of £ and the 
non-interacting semicircular DOS po as: 
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We solve for these equations using two powerful, and 
in principle exact, numerical methods: Quantum Monte 
Carlo (QMC) and Exact Diagonalization (ED) 0. A 
similar methodology was used in the study of a related 
two-band Hubbard model 14| . The QMC is a finite tem- 
perature method and is exact in the statistical sense, 
while the ED is at T = and relies on diagonalization 
of finite clusters and extrapolations to account for the 
systematic finite size effects [3|. 

We shall begin our discussion with the behavior of the 
density of states of the model in different regimes. For 
U = and at t pd — the system is an insulator since 
neither the d— orbitals at the Fermi energy can conduct 
(because they are localized), nor the p— band can con- 
duct (because it is full and well beneath the Fermi sur- 
face). At a finite hybridization t pd , however, the system 
becomes metallic, as the p and d orbitals form a par- 
tially filled band at the Fermi energy with mixed p and 
d character. In FigQ]we show the comparison of the p— 
and d— electron DOS. The one carrying most of the spec- 
tral intensity at low frequencies is the d— electron DOS 
Pd{io), since the bare atomic energy of the d— orbitals is 
at the Fermi energy. The lower panel of the figure shows 
the dramatic effect of correlations; when the interaction 
U is increased, a rather large gap opens in the DOS at 
the Fermi energy, driving the system to a Mott insulator 
state. The gap is of order U and results from the high 
energetic cost of double occupation of the d site. An in- 
teresting effect is that the size of the Mott gap Am may 
be substantially renormalized. In the inset of FigfT] we 
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FIG. 1: Density of states for the p— and d— electrons (dashed 
and solid line) for A = 1, t p d = 0.9. Upper panel shows the 
(7 = case with fi = 0.529 that gives a total occupation of 
three particles below the Fermi energy. The lower panel shows 
analytically continued QMC data for the Mott insulator state 
with U — 2, T — 1/64 and pi = 1.079 that gives n tot = 3. 
The double arrow head line indicates the large Mott gap Am- 
The inset shows Am (U) for different positions of the of the 
p— electron band e p — —6, —3, —2,-1 (top to bottom). The 
results are obtained with ED of finite clusters of N 3 sites and 
the data shown is for the limit N s — > oo. 



show the variation of Am(U) upon increasing the dis- 
tance of the p— band with respect to the d— electron en- 
ergy. Notice that the gap Am(U) is always smaller than 
the bare U, and become equal only asymptotically when 
e p — > — oo. This renormalization effect is of relevance to 
the difficult problem of the determination of the effective 
value of U in realistic ab initio calculations using DMFT 

The Mott insulator can be destabilized by either par- 
ticle or hole doping. Therefore the system has two dop- 
ing driven metal-insulator transitions. In the one band 
Hubbard model, the two transitions have the same char- 
acter, however, as we shall see, this is not the case in the 
present model. In Fig[2]we show the occupation rid and 
n p of the d and p sites for U — 2. The plateaux that 
appear around 0.7 < [i <■ 1.2 indicate the onset of the in- 
compressible Mott insulating state when correlations are 
strong. While the Mott insulator is associated with the 
energy cost of doubly occupying the local d— orbital, it is 
interesting to notice that the Mott plateau does not occur 
exactly at rid = 1 but at a higher value, which depends on 
the hybridization. The Mott state is in fact found when 
the total number of particles per unit cell is exactly equal 
to three. Thus, the object that becomes localized due to 
the strong correlations is not simply a d— electron, but 
a composite object with mixed d and p character. The 
inset of the figure shows the phase diagram in the XJ-[i 
plane, that maps the region of the Mott insulator phase 
and the transition lines to correlated metallic states. 

In FigfS] we show the DOS for the p and d electrons 
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FIG. 2: rid (solid line), n p (dashed) and ntot (dotted) as a 
function of /i, for U — 2. The data are from QMC at A = 
1, tpd = 0.9 and T = 1/64. The plateaux at n to t=2 and 
4 are band insulator (BI) states. The one at ntot=3 is the 
Mott insulator (MI). The inset shows the phase diagram in 
the U-/J, plane. The boundary lines are for T=l/20 (dotted), 
T=l/64 (thick solid) and T = (thin solid). The dashed 
thick line segment at T = 1/64 denotes the region of the 
MIT boundary where the QMC data show coexistence of a 
metal and an insulating solution. The T—0 data are from 
extrapolated ED calculations. This method is not suited for 
the study of coexistence of solutions at T = 0. The dash- 
dot line denotes the transition from a metal (M) to the band 
insulator (BI) at n to t = 4. 
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FIG. 3: Density of states of p— and d— electrons (dashed and 
solid line) for U = 2, A = 1, t pd = 0.9 and T = 1/64, as 
obtained from QMC. Upper panel: /i = 0.554, which corre- 
sponds to tiny hole doping. Lower panel: fi = 1.234, which 
corresponds to tiny particle doping. 



in the metallic states that are obtained by either particle 
or hole doping of the Mott insulator. In both cases one 
finds that the DOS clearly show the emergence of a corre- 
lated small quasiparticle peak at the Fermi energy. The 
occurrence of a narrow quasiparticle peak at the Fermi 
energy that is flanked by large Hubbard bands which are 
separated by an energy of order U is a hallmark result 
of the solution of the Hubbard model within DMFT [3| ■ 
Thus, one may be led to conclude that the MIT in the 
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PAM shares the same qualitative features. Rather sur- 
prisingly, this expectation is only fully confirmed for the 
case of particle doping, but the MIT scenario in the hole 
doped case is strikingly different. Upon particle doping of 
the Mott insulating state we have confirmed that there 
is a small region of parameters at the MIT boundary 
where two coexistent solutions, one metallic and one in- 
sulating, are found. In addition, the numerical solutions 
show critical slowing down of the convergence speed of 
the self-consistency close to the transition. These two 
features were also observed in the previous studies of the 
finite T first order MIT in the Hubbard model In 
contrast, in the hole doped case (n to t < 3) we found no 
trace of coexistent solutions down to T = 1/128. The 
solution seems to be always unique, which implies that 
the transition is of second order, i.e., through a quantum 
critical line in the U-fi plane. 

Insight on the physical reason for the qualitative dif- 
ference observed upon particle or hole doping is ob- 
tained from the behavior of the observables that mea- 
sures the d-p correlations. In Fig[J]we show the behavior 
of magnetic correlation between the d and p electrons 
(m^rrip) = ((n^f — n d[){n P i — "pi)) across the metal- 
insulator transitions. As the results show, the value of 
{m z d rrip) is sizable on the hole doped metal, while be- 
comes negligible on the particle doped side. In the parti- 
cle doped case (m^m|) is negligible because the p— band 
is already full and the extra particles mainly go to oc- 
cupy d— sites. Thus, the p— sites do not get an induced 
magnetic moment and they cannot screen the magnetic 
moment of the local d— sites. Thus the magnetic cor- 
relations develop directly among neighboring d orbitals. 
These correlations are of antiferromagnetic character due 
to the superexchange mechanism and are analogous to 
those created between neighboring sites in the Hubbard 
model case. Thus we can understand that for particle 
doping the character of the MIT in the Periodic Ander- 
son and Hubbard models is in fact the same; the p— sites 
merely allow the charge fluctuation (and thus the der- 
ealization) of the d— electrons, but they do not couple 
magnetically and do not screen the local moments. 

In striking contrast, upon hole doping the scenario is 
quite different. In this case, the system finds it is energet- 
ically favorable to create holes in the p— band and mag- 
netically bound the holes to the local moment of the cor- 
related d— site. This feature is reminiscent of the "Zhang- 
Rice" singlet formation |ll| and leads to the emergence 
of a quasiparticle peak at the Fermi energy as the singlets 
delocalize. 

In conclusion, we have investigated the doping driven 
metal-insulator transition in the Periodic Anderson 
model in the Mott-Hubbard regime. We found that the 
size of the Mott gap can be significantly renormalized by 
hybridization effects. In addition, we found that while 
both correlated metallic states at small doping show a 
small quasiparticle peak at the Fermi energy, the nature 
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FIG. 4: (m^rrip) as a function of fi, for U = 2, A = 1, t p d = 
0.9. The results are obtained from ED with a cluster of 16 
sites. 



of the MIT is qualitatively different on each side. In 
the particle doped side the quasiparticle peak is associ- 
ated with a Kondo-like resonance and the MIT shares the 
same qualitative nature of the first order transition found 
in the one band Hubbard model. In contrast, on the hole 
doped side, the quasiparticle peak is associated with the 
formation of "Zhang-Rice" singlets and the transition is 
second order. Thus, our study demonstrates that, even in 
relatively simple situations, the one band Hubbard model 
should not be automatically considered the low energy ef- 
fective model of more complicated multi-orbital systems. 
The investigation of the physical nature of the "Zhang- 
Rice" correlated metal is a very interesting problem open 
for future investigations. 

We acknowledge M. Gabay for useful discussions. 



[1 
[2 
[3, 
[4 
[5 

ir 

[§ 
[9 
[10 



N.F. Mott, Metal-Insulator Transitions (Taylor & Fran- 
cis, London, 1974). 

For a review, see M. Imada, A. Fujimori and Y. Tokura, 
Rev. Mod. Phys. 70, 1039 (1998). 

For a review, see A. Georges, G. Kotliar, W. Krauth and 

M.J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 

W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 

(1989). 

A. Georges and W. Krauth, Phys. Rev. B 48, 7167 
(1993); M.J. Rozenberg, G. Kotliar and X.Y. Zhang, 
Phys. Rev. B 49, 10181 (1994). 

M. J. Rozenberg, R. Chitra, and G. Kotliar, Phys. Rev. 

Lett. 83, 3498 (1999); G. Kotliar, E. Lange, and M. J. 

Rozenberg, Phys. Rev. Lett. 84, 5180 (2000). 

P. Limelette et al, Science 302, 89 (2003); P. Limelette 

et al, Phys. Rev. Lett. 91, 016401 (2003). 

J. Zaanen, G. A. Sawatzky and J. W. Allen, Phys. Rev. 

Lett. 55, 418 (1985). 

G. Kotliar and D. Vollhardt, Physics Today 57, 53 
(2004). 

W. F. Brinkman and T. M. Rice, Phys. Rev. B 2, 4302 
(1970). 



■5 



[11] F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 
(1988). 

[12] G. Kotliar et al., Rev. Mod. Phys. 78, 865 (2006). 
[13] A. Georges, G. Kotliar, and Q. Si, Int. J. Mod. Phys. 6, 
705 (1992). 



[14] T. Mutou, H. Takahashi and D. S. Hirashima, J. Phys. 
Soc. Jpn. 66, 2781 (1997); Y. Ono, R. Bulla and A.C. 
Hewson, Eur. Phys. J. B 19, 375 (2001). 



